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We introduce a simple, solvable model of double-barrier resonant tunneling structure which in- 
cludes the effects of electron-electron and electron-phonon scattering. The model is based on a 
generalized effective-mass equation where a nonlinear coupling is introduced to account for those 
inelastic scattering phenomena. The nonlinear term depends on one parameter which is the only 
constant in the model to be determined phenomenologically from experimental data. As an example 
of the application of the model, we discuss GaAs-Gai-^Al^As double-barrier structures. When the 
nonlinear term is considered, a sideband is observed at an energy of the order of the LO phonon 
energy in GaAs in addition to and below the main transmission peak. Fitting of the computed 
spectrum to known experimental results allows to find the value of the nonlinear coupling. As a 
consequence, we can estimate the electron-phonon coupling from the shift of the main transmission 
peak in comparison to the linear resonant tunneling. We obtain good agreement with values coming 
from more sophisticated treatments, which supports our claim that the nonlinear term consistently 
includes electron-phonon effects. Finally, we use our simple model to study magnitudes of interest 
for applications, namely I — V characteristics, which are shown to present two negative differential 
resistance peaks arising from the main transmission peak and the sideband. 

PACS numbers: 73.40. Gk, 72.10.-d, 03.40.Kf 
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I. INTRODUCTION 

Resonant tunneling (RT) through double-barrier struc- 
tures (DBS) make these systems very promising can- 
didates for a new generation of ultra-high speed elec- 
tronic devices. For instance, a GaAs-Gai-^Al^As DBS 
operating|-at THz frequencies has been reported in the 
literature III The basic reason for RT to arise in DBS is 
a quantum phenomenon whose fundamental characteris- 
tics are by now well understood: There exists a dramatic 
increase of the electron transmittivity whenever the en- 
ergy of the incident electron is close to op.e of the unoc- 
cupied quasi-bound-states inside the wello In practice, a 
bias voltage is applied to shift the energy of this quasi- 
bound-state of nonzero width so that its center matches 
the Fermi level. Consequently, the I — V characteristics 
present negative differential resistance (NDR). 

In actual samples, however, the situation is much more 
complex than this simple picture. This is so even in 
good-quality heterostructures, when scattering by dislo- 
cations or surface roughness is negligible. In particular, 
inelastic scattering is always present in real devices. Ex- 
amples of inelastic scattering events are electron-phonon 
and electron-electron interactions, in which the energy 
of the tunneling electron changes and the phase mem- 
ory is lost. The influence of these many-body effects on 
DBS has recently attracted considerable attention. Sev- 
eral approaches have been carried out to model electron- 
optical-phonon scattering in RT devices.fTTJ One thing 
that is currently clear is that electron-phonon interaction 



shifts downwards the RT peak in energy by AE — gfkJo, 
where g is the electron-phonon coupling and Tkoo is the 
LO phonon energy. It has also been noted that sideband 
transmission peaks arise, which are readily observable 
in I — V characteristics, and correspond to sequential 
rather than coherent tunneling. On the other hand, the 
effects of Coulomb repulsion in RT have also been con- 
sidered in tight-binding models of quasi-one-dimensional 
systems.! Results showed that the conductance exhibits 
several doublet structures depending upon the U cou- 
pling in the on-site Coulomb interaction between two 
electrons. Also, each peak shifts upwards to higher en- 
ergies as U increases. Therefore, we may summarize the 
available knowledge by saying that those inelastic pro- 
cesses split the main RT peak and shift its energy down- 
wards or upwards, depending on the particular scattering 
channel. 

Even with its rather satisfactory degree of success, 
many-body calculations have difficulties that, in some 
cases, may complicate the interpretation of the under- 
lying physical processes. For this reason, an alternative 
approach to inelastic scattering effects above discussed 
has recently emerged, based on nonlinear dynamics of 
single excitations in solidsH Loosely speaking, this kind 
of treatment could be regarded as similar to Hartree- 
Fock and other self-consistent techniques, which substi- 
tute many-body interactions by a nonlinear effective po- 
tential. In this paper we present one of such a nonlinear 
model, intended to account for many-body effects in RT 
through DBS within the one-particle framework, by in- 



eluding nonlinear interaction terms in the equation of 
motion. We want to make it clear from the beginning 
that we neither search for nor claim quantitative agree- 
ment between our one-particle model and existing ex- 
periments. We realize that this is a shortcut also found 
in more elaborated many-body theories and we cannot 
hope to cure this problem with a much simpler, phe- 
nomenological theory: It is obvious that if we try to de- 
scribe electron-electron and electron-phonon interactions 
within the same one-particle framework, quantitative ac- 
curacy is no longer possible. On the contrary, what we 
do claim is that our model captures the essential physics 
of inelastic effects in RT in a very simple way. The virtue 
of such an approach is that it allows to gain insight on 
the features of DBS without the burden of intensive com- 
putations providing, in an inexpensive way, a qualitative 
picture of what is to be expected in particular devices. 
This is so because our model can be easily extended to 
other systems in spite of the fact that we here present it in 
the context of a specific DBS. Therefore, our model can 
give hints leading to the development of sophisticated, 
possibly many-body calculations with much better accu- 
racy. In addition, the way we connect effective nonlin- 
earity with the physical ingredients it tries to substitute 
may be helpful for other nonlinear models in condensed 
matter physics. 

The paper is organized as follows. In Sec. II, we present 
our model, obtained by including a nonlinear coupling in 
a generalized effective-mass equation. We particularize 
it for a GaAs-Gai-^ALAs but we insist that the model 
is quite general and applicable in different contexts. The 
nonlinear coupling can be regarded as an adjustable pa- 
rameter to fit as well as possible the available experimen- 
tal results. We discuss in detail the physics underlying 
our choice for the coupling. We sketch the exact solu- 
tion of this phenomenological model, as well as the way 
to obtain the transmission coefficient as a function of 
the nonlinear couplings and the applied voltage V. For 
completeness, we include a brief discussion of the range 
of applicability of the equation and its connection with 
physical interpretations. Afterwards, Sec. Ill contains 
the main results and discussions of our analysis concern- 
ing electron transmission and I — V characteristics. We 
find the correct value of the coefficient of the nonlinear 
coupling fitting available data and, as a function of it, we 
discuss the corresponding value of the electron-phonon 
coupling. We compute the I — V characteristics of the 
device and discuss its main features, the most salient of 
them being the existence of two peaks of NDR. Finally, 
Sec. IV concludes the paper with a brief survey of the re- 
sults, and some prospects on the application of the ideas 
we have discussed. 



II. MODEL 



A. Physical grounds and definitions 

To describe our model we have chosen a typical sys- 
tem: A GaAs-Gai_ x Al x As DBS under an applied elec- 
tric field. The thickness of the whole structure is L and 
the thickness of the well is d. The barriers are assumed to 
be of the same thickness (symmetric case) but as will be 
evident below this is not a restriction of our approach. 
The structure is embedded in a highly doped material 
acting as contact, so that the electric field is applied 
only in the DBS. We focus on electron states close to 
the bandgap and thus we can neglect nonparabolicity ef- 
fects hereafter. Then the one-band effective-mass frame- 
work is completely justified to get accurate results. For 
the sake of simplicity, we will further assume that the 
electron effective-mass m* is constant through the whole 
structure. This hypothesis is related to the fact that we 
are not interested in high quantitative accuracy, although 
we note that the spatial dependence of the effective mass 
can be taken into account if necessary. 

Within this approach, the electron wave function is 
written as a sum of products of band-edge orbitals 
with slowly varying envelope- functions. Therefore the 
envelope-function ip(z) satisfies a generalized effective- 
mass equation (we use units such that energies are mea- 
sured in effective Rydberg (Ry*) and lengths in effec- 
tive Bohr radius (a*), being IRy* = 5.5 me V and I a* = 
100 A in GaAs) given by 

- iM*) + [V(z) X b(z) - eFz] %l>(z) = E i>(z), (1) 

where V(z) is the potential term which we discuss be- 
low, F is the electric field, and Xb(z) is the characteristic 
function of the barriers, 



if < z < (L - d)/2, 
X b(z)= { I, if (L + d)/2<z<L, 
otherwise. 



(2) 



We now specify our model by choosing what is the po- 
tential term V(z). In order to do that, let us first consider 
the physics we are trying to represent with this term. The 
DBS can be regarded as an effective medium which re- 
acts to the presence of the tunneling electron, leading to 
a feedback mechanism by which inelastic scattering pro- 
cesses change the RT characteristics of the device. It thus 
follows that V{z) must contain nonlinear terms if it is to 
summarize the medium reaction which comes from the 
electron-electron and electron-phonon interaction. The 
simplest candidate to contain this feedback process is the 
charge density of the electron, which is proportional to 
| , i/'(^)| 2 - In our model, we neglect higher order contribu- 
tions and postulate that the potential in Eq. (fil) has the 
form 



V{z) = V b [l + &\^{z){< 



(3) 



where VJ, is the conduction band-offset at the interfaces, 
and all the nonlinear physics is contained in the coeffi- 
cient a which we discuss below. 



There are two factors that configure the medium re- 
sponse to the tunneling electron. First, it goes without 
saying that there are repulsive electron-electron Coulomb 
interactions, which should enter the effective potential 
with a positive term proportional to the charge, i.e., the 
energy is increased by local charge accumulations, lead- 
ing to a positive sign for a. On the other hand, in polar 
semiconductors, the electron polarizes the surrounding 
medium creating a local, positive charge density. Hence 
the electron reacts to this polarization and experiences 
an attractive potential (which implies a negative). This 
happens, for instance, in the polaron problem in the weak 
coupling limit, which becomes valid in most semiconduc- 
tors, and where it can be seen that the lowest band en- 
ergy state decreases^ It is then clear that in principle 
any sign would be equally possible for this coefficient if 
a is to represent the combined action of the polarization 
of the lattice along with repulsive electron-electron inter- 
actions. Intuitively, however, it is most realistic to think 
that a will be negative, because a positive nonlinear in- 
teraction would arise from negative charge accumulation 
in the barriers, which is not likely to occur. In fact, we 
discuss below mathematical reasons imposing that a has 
to be negative as expected, allowing for it to be positive 
only if it is small and the barriers themselves are very 
narrow. It can be argued following this line of reason- 
ing that charge accumulation is expected inside the well, 
and that we should include a nonlinear term in the well 
to account for that. We have studied this possibility as 
well and we found that such a term is irrelevant, as we 
explain in detail in Sec. III. 



B. Analytical results 

In the preceding subsection, we have defined our 
model, and the only remaining thing is to fix the value 
of a. This we address in the next section. We now work 
starting from Eq. (fil) with the definition in Eq. (0) to 
cast the equations in a more tractable form. For sim- 
plicity, and because we are interested in intrinsic DBS 
features, we consider that the contacts in which the struc- 
ture is embedded behave linearly. Therefore, the solution 
of Eq. (||) is a linear combination of traveling waves. As 
usual in scattering problems, we assume an electron in- 
cident from the left and define the reflection r and trans- 
mission t amplitudes by the relationships 
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z<0, 

z > L, 



(4) 



where fc 2 



E,kl 



E + eFL, and A is the incident wave 



amplitude. Now we define ?p(z) — A<fi(z), and a = <S|AP. 
Notice that a is a dimensionless parameter. Using Eq. (|l|) 
we get 



<t>zz{z) + {V bXb (z) [1 + «|0(z)| 2 ] - eFz - E) cf>(z) = 0, 

(5) 



To solve the scattering problem in the resonant struc- 
ture we develop a similar approach to that given in 
Rcf. FLO. Since (f>(z) is a complex function, we take 
4>{z) = q(z) exp[z7(z)], where q(z) and "f(z) are real func- 
tions. Inserting this factorization in Eq. (0) we have 
j z (z) = q~ 2 {z) and 

- qzz(z) + ^-r + [V b Xb(z) - eFz - E] q(z) + 
q A {z) 

+aV bX b{z)q 3 (z) = 0. (6) 

This nonlinear differential equation must be supple- 
mented by appropriate boundary conditions. However, 
using Eq. (0) this problem can be converted into a ini- 
tial conditions equation. In fact, it is straightforward to 
prove that 

q(L) = fc- 1/2 , q z (L) = 0, (7) 

and that the transmission coefficient is given by 
4fc 9 2 (0) 



l + 2M 2 (0) + fc V(0) + g 2 (0)<z, 2 (0) 



(8) 



Hence, we can integrate numerically (0) with initial con- 
ditions (ffl) backwards, from z = L up to z — 0, to obtain 
q(0) and q z (0), thus computing the transmission coeffi- 
cient for given nonlinear coupling a, incoming energy E 
and applied voltage V — FL. 

Once the transmission coefficient has been computed, 
and recalling that contacts are linear media, the tunnel- 
ing current density at a given temperature T for the DBS 
can be calculated within the stationary-state model from 



m = 



m*ekBT 
2n 2 h 3 



T(E,V)N(E,V)dE, (9a) 



where N (E, V) accounts for the occupation of states to 
both sides of the device, according to the Fermi distribu- 
tion function, and it is given by 

N(E,V) = lJ 1 + ^-WI V (9b) 



I + cxp[(£ F -E- eV)/k B T 
where ks is the Boltzmann constant. 

C. Sign of the nonlinear term 



At this point, we are in a position to give mathemat- 
ical reasons why a (or a) should be negative by using 
Eq. (P). This equation has the form 

- qzz(z) + -a^r + fl(z)q(z) + f 2 (z)q 3 (z) = 0, (10) 

where fi(z) and J2{z) are well-behaved functions. Let us 
now consider the discrete version of this equation, with 



the second derivative discretized in the usual way; denot- 
ing q n = q(z = nAz) and f in = fJnAz), i = 1,2, with 
Az being the integration step, Eq. (10) can be rewritten 
as 



Qn+i — 2<7„ 



(Az) 2 



~o + flnQn + hnq n 
9n 



(11) 



Notice that in this expression the sign of fin is the 
same as the sign of a and a. Let us now consider Eq. ( p"l| ) 
for large q n ; this is general, because if q is small the term 
q~ 3 will make it grow quite quickly. In this limit, Eq. (p/fl) 
can be approximately replaced by 



Ag„_! = Aq n + {Azf f 2n ql, 



(12) 



where Aq n = q n — q n +i, and we have cast the equation in 
this fashion because it is to be integrated backwards. Re- 
calling the initial conditions (M), if iV is the total number 

of grid points, we have gjv = 5/v-i = k L 
fore, we see that Agjv-i = 0, and 



1/2 



> 0. Thcrc- 



Aq N _ 2 = (Az) 2 f 2n k 



-3/2 



(13) 



We thus see that if f 2n is positive (and hence a) the first 
increment is positive, and so are the subsequent ones, 
leading to a exponential divergence of q, whereas if f 2n 
and a are negative, the increment is negative, and q de- 
creases until the q~ 3 starts being relevant again. In this 
last case, it is possible that q reaches an equilibrium due 
to the balance of the two cubic terms, which was not for 
a > 0. This is in fact seen in the numerical integration 
of Eq. (@), where q rapidly diverges if a is positive, un- 
less, of course, a is positive but very small and/or the 
barrier (the region for a to influence q) is very narrow. 
In this last situation, however, the effect of a becomes 
negligible. 

Physically, this can be understood as follows. Elec- 
trons impinge on the barrier from outside, and begin to 
tunnel through it, their wavefunction being real and ex- 
ponentially increasing or decreasing. If a is positive, then 
if there were any charge density in the BDS, this would 
become even more repulsive, and the wavefunction would 
diverge even faster (numerically one always see the expo- 
nentially growing part, of course). This instability is not 
present in the opposite case, where a negative a helps 
the electron tunnel across the barrier. 



III. RESULTS AND DISCUSSIONS 
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In our calculations we have considered a double-barrier 
GaAs-Gao.65Alo.35As structure with L = 3d = 150 A. 
The conduction-band offset is Vb — 250 meV. In the ab- 
sence of applied electric field and nonlinearities, there 
exist a single, very narrow resonance r ~ 1 below the 
top of the barrier, with an energy of 80.7 meV. 



FIG. 1. Transmission coefficient r as a function of the elec- 
tron energy for a = —7.6 x 10~ 3 , for different values of the 
applied voltage (a) volts, (b) 0.029 volts, and (c) 0.105 volts. 
For comparison, dashed lines indicate the results for a — 0. 



Hence the well supports a single quasi-bound state. 
When there exists an applied voltage, the energy of 
the quasi-bound state level is lowered and a strong en- 
hancement of the current arises whenever the Fermi level 
matches this resonance, thus leading to the well-known 
RT phenomenon. This picture changes dramatically if a 
small amount of nonlinearity appears in the structure. 

Figure |l] shows the transmission coefficient as a func- 
tion of the incoming energy for different values of the ap- 
plied bias, when a = —7.6 x 10~ 3 . This value is the one 
that best fits the experimental splitting that we discuss 
now. At zero bias, shown in Fig. 0(a), there exist two well 
differentiated peaks, the main one centered at 78.2 meV 
and the sideband centered at 41.5 meV. The main peak 
not only shifts to lower energies but also broadens, in 
a similar fashion to what happens in RT when inelastic 
effects arise. The energy separation between the main 
peak and the sideband is 36.7 meV, very close the energy 
of LO phonons in GaAs (hu>o — 36.3 meV), which was our 
requirement to choose a. Therefore, the self-interaction 
we introduce in our model gives results similar to those 
found in dealing with the full electron-phonon Hamil- 
tonian in DBS, a where the sideband is due to phonon 
effects at inelastic channels. Interestingly, at the side- 
band the transmission probability is close to 0.2, which 
agrees with more sophisticated treatments.u This agree- 
ment reinforces our hope to be treating in a proper way 
these inelastic scattering effects. Furthermore, to check 
the validity of our model, we can estimate the electron- 
phonon coupling from the downwards shift of the main 
peak when nonlinearities are introduced. This shift is 
AE = 2.8 meV, and using the value of Hujo = 36.7 meV 
obtained previously we get g = 0.076, in rather good 
agreement with the value g = 0.10 proposed by Cai et all3 
We note, however, that this result is not far from g = 0.03 
given by Wingreen et alB In fact, both many-body re- 
sults enclose ours, and our calculation is not clearly fa- 
vorable to any of them. 

We now consider the effect of a bias imposed on the 
DBS. When there exists an applied voltage, the trans- 
mission peaks also shifts to lower energies due to the 
presence of the linear potential, as shown in Fig. 0(b). 
This behavior is similar to that found in linear (coherent) 
RT. For instance, these peaks are located at 62.9 meV 
and 15.2 meV for V — 0.029 volts. A further increase of 
the applied voltage makes the sideband to disappear [see 
Fig. |l](c)]. The reasons for that will be explained below. 

In order to gain insight on the nonlinear RT in DBS, 
we can rewrite Eq. (ph as follows 



4>zz{z) + V eff (z, E) cf>(z) = E <f>(z), 



(14a) 



where we have defined an effective potential as follows 

V eff (z,E) = V bXb (z)[l + aq 2 (z)] - eFz. (14b) 

Thus (jlj) is a Schrodinger-like equation for an effective 
potential due to nonlinearity plus the linear potential and 
the built-in potential of the DBS. This effective potential 



depends not only on z but also on the incoming energy 
E through the function q(z). Since the envelope func- 
tion changes under RT conditions, and also q(z) accord- 
ingly, it should be expected that V e ff(z,E) undergoes 
severe variations whenever E is close to one of the RT 
peaks. This is indeed the case, as shown in Fig. |[ where 
V e ff(z, E) is displayed at zero bias. Notice that nonlinear 
effects have negligible effects on the shape of the effec- 
tive potential in the right barrier, besides a slight band 
bending at the interface z = (L + d)/2 at low energies. 
However, the potential in the left barrier region differs 
significantly from the original square-barrier shape. Out 
of resonances, the effective barrier height at z = is lower 
than Vb, whereas at resonances it takes the value ~ V&. 
Hence the effective potential presents two local maxima 
in the plane z = as a function of the incoming energy 
E, matching the values of the main resonance and the 
sideband above discussed. From Fig. @ it is clear that at 
the energy of the main resonance the effective potential is 
quite similar to the built-in DBS potential, just produc- 
ing a small shift of the quasi-bound-state in comparison 
to the linear RT process. Concerning the sideband, the 
effective potential presents a deep minimum at the in- 
terface z = (L — d)/2, which causes the lowering of the 
quasi-bound-state, thus explaining the origin of this lower 
RT peak. Indeed, the fact that this added well is respon- 
sible for the peak is confirmed by the observation that the 
disappearance of the peak as the field increases coincides 
with the vanishing of the well due to the bias; before 
that, the resonance responsible for the peak is moving 
towards E = where it can be seen that the depth of the 
secondary well vanishes. We thus provide a complete, 
coherent picture of the tunneling phenomenology. 
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FIG. 2. Effective potential V e jf as a function of z and the 
incoming electron energy E at zero bias, for a — —7.6 x 10~ 3 . 

We now have to comment the j — V characteristics, 
computed from (0). We have set two different tempera- 
tures (77 K and room temperature) and compared these 
curves with those obtained in linear RT. The Fermi en- 
ergy was Ep — 27.5 meV. Results are shown in Fig. ||. 



In all cases we obtained clear NDR signatures. The lin- 
ear RT shows a single NDR peak whereas self-interaction 
causes the occurrence of a second peak at a lower voltage, 
clearly related to the sideband in the transmission coef- 
ficient. On increasing temperature both peaks merges 
into a single one because of the broadening of the Fermi- 
Dirac distribution function. It is important to notice 
that the j — V characteristics collect all the main features 
found in the transmission coefficient, namely shift down- 
wards of the main RT peak, its corresponding broadening 
due to inelastic effects and the appearance of the side- 
band, shown in Fig. [| as a smaller NDR signature at 
V = 0.05 volts. 

Finally, some words are in order about possible non- 
linearities in the well, as we announced in Sec. II. In 
principle, charge may accumulate in the well depending 
upon doping degree, and that is expected to induce a self- 
repulsive interaction of the electronic cloud, overcoming 
the attractive part due to polarization of the medium. So, 
if we are considering that nonlinearities are zero in the 
well we could argue that one effect compensates the other 
and then there is no nonlinear term to be added. Admit- 
tedly, that would be a very rare phenomenon, because 
exact compensation of both forces (or any interaction in 
nature) would be serendipitous. The question then arises 
as to whether our model is a good one in the sense that 
it is robust against new terms accounting for inelastic 
scattering in the well. We checked that by introducing 
another term in the potential, given by V'(z) = /3|?/>| 2 
for (L - 2)/2 < z < (L + d)/2. The calculations can 
be carried out much in the same way we have discussed, 
and the results obtained by numerical integration of an 
equation very close to Eq. (|j). We studied a range of 
values for j3 = (3\A\ 2 of the same order as those of a: (3 
positive would indicate that Coulomb effects dominate, 
whereas f3 negative would be the same as in the barriers, 
a representation of the effects of polarization. The results 
were very much satisfactory, because the features of the 
model did not change, and only minor quantitative mod- 
ifications were found. We can then conclude that leaving 
the possible nonlinearity of the well out of our model is 
not a limitation and that it is not necessary to take it 
into account. 



3:M) 





0.00 



(a) 






J A 
{ f 


- 


// 




J \ 


j 




/ / 

J L -v 


'•' 





0.05 0.10 0.15 

Applied voltage (Volts) 



50 - 



0.20 




0.00 0.05 0.10 0.15 

Applied voltage (Volts) 



FIG. 3. Computed j — V characteristics for Ef — 55meV, 
a — -7.6 x 10~ 3 at (a) T — 77 K and (b) room temperature. 
For comparison, dashed lines indicate the results for a — 0. 



IV. CONCLUSIONS 

In conclusion, we have presented a nonlinear effective- 
mass equation to take into account in a simple way 
electron-electron and electron-phonon interaction in RT 
through DBS. Our model represents those effects by 
including a self-coupling of the wavefunction with the 
charge cloud inside the barriers, whose strength is the 
only adjustable parameter. By comparing the results of 
the model to experimental data on the distance between 
the main transmission peak and the sideband we fix the 
value of that constant. 



With the same value we found transmission probabil- 
ities and electron-phonon coupling values very close to 
those arising from many-body calculations, which indi- 
cates that our model consistently reproduces their fea- 
tures. We have also been able to discard any possible 
nonlinearity in the well, which shows that our under- 
standing of the physics of this problem is in the right 
direction. Wc have to stress that the low values obtained 
for the nonlinear coupling, a = —7.6 x 10 -3 , do agree 
with our purpose of offering an alternative description: 
Had we obtained large values of a, it would be very dif- 
ficult to understand how those magnitudes could arise 
from what in principle are perturbative effects of the lin- 
ear description. In addition, large values of a would have 



immediately posed the question about the role that soli- 
ton formation and transmission could be playing along 
with RT in the DBS (see Refs. || and [n] as well as 
references therein). Therefore, our model satisfactorily 
achieves the goals we had in mind, a great success if its 
simplicity is taken into account. 

Further extensions of the present work to study non- 
linear dynamical response of DBS on external ac bias 
would be of great interest to shed light on plated prob- 
lems like bistability,E3 npise characteristics ]13 and RT at 
far infra-red frequenciesE.il under the influence of inelas- 
tic scattering channels as those described here. Besides 
that, the above mentioned highly nonlinear limit could 
also be interesting, for nonlinearity is always susceptible 
to give rise to new and unexpected features. If these new 
features were seen in our model it would be a very excit- 
ing development. If materials with suitable characteris- 
tics were found, and that is to be expected, experiments 
could be made to check the predictions: If the model were 
wrong, that would establish its range of validity, whereas 
if the predictions were correct, this work could pave the 
way to a new family of devices and applications. 
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